Numerical study of current distributions 
in finite planar Hall samples with disorder 
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A numerical study is made of current distributions in finite- 
width Hall bars with disorder and some theoretical observa- 
tions are verified. The equilibrium current and the Hall cur- 
rent are substantially different in distribution. It is observed 
that in the Hall-plateau regime the Hall current tends to con- 
centrate near the sample edges while it diminishes on average 
in the sample interior as a consequence of localization. The 
edge states themselves scarcely affect both the equilibrium 
and Hall currents in the plateau regime. The sample edge 
combines with disorder to efficiently delocalize electrons near 
the edge while the Hall field competes with disorder to delo- 
calize electrons in the sample interior. A possible mechanism 
is suggested for the breakdown of the quantum Hall effect. 

73.40.Hm,73.20.Jc 



I. INTRODUCTION 

The current distribution in a Hall sample is a subject 
of interest-|De£taining to the foundation of the quantum 
Hall effecrali-3 (QHE). The standard explanations of the 
QHE are based on a picture that the Hall current is car- 
ried by a finite fraction of electron states in the sample 
bulk that remains extended in-the presence of disorder; 
there localization is essentialffd for the formation of vis- 
ible conductance plateaus. An alternative explanation is 
based on a pictureolij in which the electron edge statescl 
are taken to be principal carriers of current, and leads to 
the concept of edge channels .-which turns out useful in 
summarizing observed results i2_3 Still it is not quite clear 
experimentally whether the curreatrflpws only along the 
sample edges; some experiments Jij~li3 though indirectly, 
suggest edge-channel widths much larger than the width 
naively expected from the edge-state picture of the QHE. 
It is important to clarify Japp. these apparently different 
pictures could be relatedJlJci Theoretical consideration 
on current or potential distributions has so far focused 
on clean samplescj rather than disordered samples E3 

The purpose of this paper is to present a numerical 
study of current distributions in small Hall bars with dis- 
order and to verify some theoretical observations. We 
shall handle samples which accommodate as many as 
280 electrons; this makes it possible to disclose features 
of the current-distributions not observed in earlier nu- 
merical worker on smaller samples. In particular, we 
shall find evidence that in the Hall-plateau regime the 



Hall current tends to concentrate near the sample edges 
while it diminishes on average in the sample interior as a 
consequence of localization. The edge states themselves 
scarcely affect the distributions of both the equilibrium 
and Hall currents in the plateau regime. 

We shall use a uniform Hall field to detect the current- 
carrying characteristics of Hall electrons. One certainly 
has to treat the problem of an electrostatic potential self- 
consistently, in which chargeS|-adll be redistributed as 
pointed out by Chklovskii et al£3 What we propose here, 
however, is that, even before the self-consistent potential 
is used, there exists an important physics of redistribu- 
tion of the Hall current via disorder, and we believe that 
even after the inclusion of the selfconsistency the essen- 
tial feature will remain because the key mechanism is 
competition of the electrostatic potential with disorder 
rather than its distribution, as we shall see. 

In Sec. II we present a theoretical framework for our 
numerical simulations. In Sec. Ill we discuss the influence 
of isolated impurities on electron states. We handle a 
Hall sample with an array of impurities in Sec. IV and 
samples with random impurities in Sec. V. Section VI is 
devoted to a summary and discussion. 



II. FORMALISM 



Consider electrons confined to a strip of length L x and 
width L y (or formally, a strip bent into a loop of circum- 
ference L x ), described by the Hamiltonian: 



H = H a + U(x,y)~eA Q (y), 

H = ^{£V y + (i/^)(v-yo) 2 }, 



(2.1) 
(2.2) 



written in terms of u> = eB /m, the magnetic length £ = 
1 / Vei? and ?/o = p x /(eB). Here the Landau-gauge vector 
potential (—By,0) has been used to supply a uniform 
magnetic field B normal to the plane. We take explicit 
account of the edge y = 0, where the wave function is 
bound to vanish; the prescription for the other edge y = 
L y will be specified later. 

We shall detect the current j x (x, y) flowing in the pres- 
ence of a Hall potential Ao(y) and an impurity potential 
of the form 



U(x,y) = 22\i 8(x - Xi)5(y - yi), 



(2.3) 
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which consists of short-range impurities of strength \ 
distributed over the sample at position (xi, yi). For sim- 
plicity, we here take Ao(y) — —yE y which produces a 
uniform Hall field E y in the y direction. 

The kinetic term Ho describes the cyclotron motion 
of an electron, which gives rise to Landau levels \N) = 
\n, yo) labeled by integers n = 0, 1, 2, • • •, and yo = £ 2 Px, 
with wave functions of the form 



(x,y\N) = (L x ) 



-1/2 Jx. 



p "0n(y;2/o)- 



(2.4) 



The dia.(y; yo) are givenEl by the parabolic cylinder func- 
tionscj D v (\[2(y — yo)/£)i which for electrons residing 
in the sample bulk are reduced to the usual harmonic- 
oscillator wave functions. Setting <f) n — at y = fixes 
the energy eigenvalues 



e N = u{v n (y Q ) + i} 



(2.5) 



as a function of yo for each n. The spectra v n (yo) rise 
sharply only for yo ^ 0(£) and recover integer values n 
f° r Vo ^ 0{& few £) p for the n = level, voiuo) decreases 
from 1 to 0.003 as yo varies from to 2.5£. 

The presence of the sample edge y = has effec- 
tively generated strong potentials ujv n (yo) that confine 
electrons to a region yo > in yo space and that 
drive electrons along the edge. Indeed, near the sample 
edge an electron state (n,yo) travels with velocity v x « 
Lv£ 2 v n '(yo) ~ — much larger than the field-induced, 
drift velocity E y /B. Classically these edge statesfflliMj 
are visualized as electrons hopping along the sample edge 
and carry current even in equilibrium. 

The 4> n (y;yo) are highly localized around y ~ y cm = 
yo - £ 2 v n (yo) with spread Ay ~ 0(1); see Eq. ( |2~9| ) be- 
low. With (x,y\N) taken to be periodic in x, yo = £ 2 p x 
is labeled by integers k: 



y Q = (2TT£ 2 /L x )k = yo[k}. 



(2.6) 



We shall henceforth use this discrete label k rather 
than yo to refer directly to each electron state, denote 
\N) = \n,yo[k]) simply as \n, k) and use the normaliza- 
tion (jv|jv'> - <W<W- 

For numerical calculations it is advantageous to handle 
the Hamiltonian Hnn' = (iV|H|JV') in N = (n,k) space 
rather than in real (cc, y) space: 

H NN , = uj{v n (yo) + \}5nn'+ U NN >+ eE y y NN >, (2.7) 
with 

£W = V A ie - l ^-^)^/ £2 0„(y i; yo)0„' (y l ;y'o), (2-8) 

% 

VNN' = {Vo <W + £ Ynn' (Vo)} 8kk> , (2.9) 

where U NN , = (N\U\N'), y NN , = (N\y\N'),y = y [k] 
and y' = yo[k'A- The matrix elements Y nn >(yo) are com- 
pletely knownEil through the spectra v n {yo)- 



with v n ' = d Vo v n (yo)\ in particular, Y nn (y ) = 
—£vn(yo)- For yo S> 0(£), Y nn > are reduced to constant 

y/n'/2 6 n +l,n' 



matrices Y^ n T"' = ^n'/2 5 n+ i^ + \fn/2 5 n , nl+l . 

When the disorder Unn' is weak compared with the 
level gap, i.e., \Xi/2Tr£ 2 \ < u>, one can diagonalize Hnn' 
with respect to Landau-level labels (n, n') by a suit- 
able unitary transformation H — \VHW~ 1 so that 
the diagonal component (N\H W \N') = 8 nn > H™(k,k') 
describes each impurity-broadened Landau subband. In 



particular, with a simple 0(U) choice W 



with 



iA 



NN' 



eE y yNN' + Unn 1 



(for n ^ n'), (2.11) 



u(v n [k) - v n t[k']) 
the subband Hamiltonian H^ l (k, k') reads 

HZ(k,k') ={w(v n [k]+ \) + eE v {yo-£ 2 v n '(yo)}}5 k k> 
+U mL (k,k') + 0(U 2 /Lo), (2.12) 

where we have set v n (yoW\) — * v n [k] and Unn 1 —> 
Unn'(k,k'). 

Some care is needed here. The choice ( [2.11 ) becomes 
inadequate for v n — v n * ~ 0; such degeneracy among 
different subbands is inherent to edge states. Fortunately, 
to the order we work (i.e., QLU)), there is a practical way 
out of such complications :cj In what follows we simply 
focus on the n = subband and consider only electron 
states with yo = {2it/L x )k > (i.e., k — 1,2, •••) that 
are free from intersubband degeneracy. 

The Hamiltonian governing the n — subband of our 
interest is H Q ^(k,k') with fc, k' — 1, 2, • • • , N max , where 
JV max is the maximum number of electrons to be retained 
in numerical simulations. Associated with each eigen- 
state \a) of Hqq is an (7V max -component) wave function 
ak = (0, A;|ai), which is found as an eigenvector by diag- 
onalizing H^(k,k') numerically. 

The x-space wave functions (x|W _1 |a) are readily re- 
constructed from these a k - 



(x\W- 



n . k 



[x.\n, k) (Son a k + 0%) 



where (x| = (x,y\ for short and 

_ _1 e£E y Y n o(yo)S k k' + U n o(k 7 k') 
Pk ~ u v n [k]-vo[k'] 



(2.13) 



aw (2.14) 



for n > 1 and /3° = 0. The spatial distribution of the 
mode a is given by 

p^(x,y) = \(x\W- 1 \a)\ 2 (2.15) 
and the current density j x = euj ^(y — | £ 2 V x )^ by 
jt\x,y) = ecuRe [(a\W\x) (x\(y - yo)^" 1 !")] , (2.16) 



where (x|(y — yo)!^ 1 |a) is given by Eq. (2.13) with 
Y nm = (-l) n+m £\v n 'v m '\ l l 2 /[l - (v n - v m f\ (2.10) (x|n,fc) replaced by (y - y [k])(x\n,k). To extract the 
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Hall-current component out oi-q 
vary E y by a small amounted 
change 



(a) 



numerically, we shall 



SE„ and calculate the 



V) = &\x, y\E y + SE y ) - i a \x, y\E y ). (2.17) 

For the simulations to be given in Sees. Ill and 
IV we have constructed the normalized wave functions 
0n(y;yo) = C n (y )D v (->/2(y - yo)/l) with v -> v n {y Q ) 
as real functions of y and yo by determining the spec- 
tra ^ n (j/o) an d overall normalization C n (yo) numerically 
to accuracy of 10 digits for y > and n =0,1, and 2. 
For each given configurati on of impurities we shall diag- 
onalize H$(k,k') of E q. (|2 12|) to 0(C/), construct the 
x-space wave functions (2.13) through the eigenvector afc, 



and calcul ate th e current ji ' of Eq. ( p. 16 ). For the wave 



functions ( 2.1 3| ) we shall take only the n = <-* n = 1 
transitions (/3F) into account; actually we have in some 
cases included the n = <-> n = 2 transitions as well and 
found no appreciable change in current distributions. 



III. ISOLATED IMPURITIES 

In the absence of impurities all electron states \n,yo) 
are plane waves extended in x, though localized in y with 
spread Ay ~ 0{£). Such spatial characteristics of elec- 
tron states are modified in the presence of disorder, and 
electrons tend to be confined in finite domains of space 
(or trapped by impurities). In this section we examine 
the influence of isolated impurities on electron states. 

For most of the cases to be discussed in this paper we 
set e£E y /uj = — 1/10 5 , which corresponds to a Hall field 
of strength E y - 0.1 V/cm under a typical setting £ ~ 
100 A and uj ~ 10 meV. [The choice E y < is taken sim- 
ply for convenience so that the y — edge is the "upper" 



edge. We have examined the E y > cases as well and 
found no essential change in our conclusion.] To calculate 



JHaii 01 E q- we set SE v/ E v = V 100 and establish 



one-to-one correspondence between the eigenfunctions of 
the SE y = and SE y ^ cases by direct examination 
of their overlap. In what follows we shall employ dimen- 
sionless impurity strengths defined by \i/2ir£ 2 — s, u>, 
and measure energy relative to \ui. 

Consider a Hall sample of length L x = \/2tt£x 10 » 25£ 
and width L y = y/2n£x9 « 22. 5£, supporting 90 electron 
states in the lowest level n = with < yo < L y , or 
yo = (2n£ 2 / L x )k with k = 1, 2, • • • , 90. The sample edge 
y = is a sharp edge while the region y ~ L y is left 
impurity-free. 

Let us place three impurities of strength (s 1; s 2 , s 3 ) = 
(0.1,0.1,-0.1) at y = (2.5£, 10£, 18£) along the line x = 
20£ on the sample and two more impurities of strength 
(si,s 5 ) = (0.05,-0.05) at y = (8£,17£) along x = Wl 
These impurities are sufficiently far apart from each other 
and can practically be regarded as isolated impurities. 
Four of them capture electrons, leading to four localized 



states of energy e; w s,u with i = 2, 3, 4, 5, as seen from 
the density distribution p^(x,y) in Fig (a). 

To get a definite idea of localized states it is useful 
to refer to the case of an isolated impurity of strength 
s = X/(2tt£ 2 lj) placed at x — y = on an infinite plane 
with E y = 0. In the n — level one localized state of 
energy e«sw arises with a wave function of the form 



loc 



1 



V2lr£2 



r (r 2 —2ixy) 



1 - s(l 



2P-' 



(3.1) 



where r 2 — x 2 + y 2 ; the rest of the n — states have 
no energy shift, e = 0. Th e O(s ) correction comes from 



the 01 cx U w term of Eq. (fulj) . Numerically Eq. Q 



accounts for the four localized states in Fig. |l| (a) very 
well. A closer look at Fig. [l] (a) reveals that a localized 
state gets somewhat narrower for an attractive impurity 
(s < 0); thus attractive impurities are more efficient in 
trapping electrons. 

A Hall field causes a drift of an orbiting electron and 
works to delocalize it. A simple estimate of energy cost 
suggests that an isolated impurity of strength s would 
fail to capture an electron for 



i\w= \\\/{2n£ 2 ) < e£\E y 



(3.2) 



Indeed, we have verified that an impurity as weak as s = 
10~ 4 can capture an electron while it fails for s = 10~ 5 . 

Similarly, a confining potential u; 1/0(2/0) works to delo- 
calize electrons. The edge states driven by an "effective" 
field as strong as e£Ey ~ <-o£uo ~ w are robust against 
disorder and remain extended.El Indeed, Fig. [I] (a) demon- 
strates explicitly that an impurity of strength s\ = 0.1, 
situated (at y — 2.5£) close to the edge, fails to trap any 
electron. (Note that £v ' ~ -0.011 for y Q = 2.51) 

Figure [l] (b) shows the current distribution jx a \%, y), 
plotted in units of — euj/L x , for each of the five states. 
There each localized state is accompanied by a circu- 
lating current distribution that reflects the underlying 
counterclockwise cyclotron motion of an electron. These 
localized current distributions are not very sensitive to 
impurity strengths Sj. In contrast, as seen from Fig. |l| 
(c), the Hall current density jj?)n(x, y) localized about 
an impurity exhibits characteristic local variations, which 
depend on the nature, attractive or repulsive, of the im- 
purity and which, somewhat unexpectedly, get amplified 
in magnitude and range as the impurity becomes weaker. 

These features of ji^ rand are readily understood 
from the current densityE3 about the localized state (JO): 



4tt£ 2 



pE 

e -f^L{y 2 -£ 2 ) 
s UJ 



(3.3) 



where, for simplicity, the 0{E y ) correction has been cal- 
culated by treating the Hall potential as a perturbation 
solely to the O(s ) part of -!/> loc . Note that the net current 
vanishes, J dyj 1 ° c (x,y) = 0. 

For a filled subband the effect of impurities disap- 
pears from the Hall-current distribution jnaa(x,y) = 
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12a 3na\\{ x i V)> which then recovers the same distribution, 
uniform in x, as in t he impurity-free case; this is verified 
by use of Eq. ( 2.13 ). In view of this fact, a vacant lo- 
calized state about a repulsive (A > 0) impurity and an 
occupied localized state about an attractive (A < 0) im- 
purity act alike, as far as jHaii(^,y) is concerned. For 
such a state jnaxi(x,y) flows on both sides of an impu- 
rity in the same direction as if it avoids the impurity 
center [like the j/) about an attractive impurity 

in Fig 0(c)]; here we already see the potential that each 
impurity expels the Hall current out of its center. 



IV. AN ARRAY OF IMPURITIES 



Let us next consider a Hall sample of length L x 



2n£ x 12 m 3(M and width L v 



2n£ x 9 w 22. U, 



on which 63 impurities of equal strength s = 0.2 are 
placed at "even" points (a^y,-) = (2(i + 2)£,2(j + 1)1) 
with i = 1, 2, • • • , 9 and j = 1, 2, • • • , 7 within the domain 
6£ < x < 2U and U < y < 16£. We have distributed 
impurities in a regular pattern so as to observe the coop- 
erative effects of many impurities in an amplified fashion. 
The remaining region 16£ < y < L y is left impurity-free 
so that it can simulate a gentle edge region supporting 
extended electron states. 

This sample supports 108 electron states in the low- 
est subband n — with < yo < L y . We arrange 
them in the order of descending energy and use the num- 
ber of vacant states, N v = 108 — N e , to specify the fill- 
ing of the subband. For convenience we start with the 
filled subband N v = and study the currents this sub- 
band supports, j x (x, y) = J2 a 3x ( x > v) and jHallfo v) = 
S« iHail( x ' y)i ^A removing electrons one by one. In 
Fig. g (a) we plot the total currents J x = J dyj x (x,y) 
and Jnaii = J dy jnedi{x, y) as functions of N v , in units of 
0.05 x (—euf£/L x ) and — (eSE y / B)(l/L x ), respectively. 

Direct inspection of the spatial distribution p^ a \x,y) 
for each state reveals that most of the first 48 states 
are localized, except for the edge states of energy « 
0.78 uj,e 2 « 0.59 w,e 3 « 0.44w,e 2 2 « 0.31 u, etc. The 
edge states are readily identified by a large amount of 
current they carry. Actually they are responsible for 
remarkable steps in J x , which emerge each time a new 
edge state becomes vacant with increasing N v . In con- 
trast, they leave little effect on Jnaii in Fig. || (a), where 
a Hall plateau persists for N v < 48. In fact, the first 
three edge states combine to carry a larger amount of 
current j x than the rest of states combined, but the 
amount of Hall current they carry is less than a single 
unit (-e5E/B)(l/L x ). 

This leads to an important observation: In the regime 
supporting clear Hall plateaus, the fast-traveling edge 
states are practically vacant and scarcely contribute to 
the currents j x {x,y) and jH_aii(x,y). The edge states, 
traveling in the negative x direction, give rise to a promi- 
nent edge current (of paramagnetic nature) for the filled 



subband (N v = 0). However, as soon as the first three 
edge states become vacant, the orbital diamagnetic cur- 
rent takes over and j x (x,y) reverses its direction in the 
edge region, as seen from the current profile across the 
sample width at x = 16 £ in Fig. || (b). It is now the 
electron bulk states that are responsible for the promi- 
nent current distributions flowing in opposite directions 
at opposite edges for 4 < N v < 50. 

As N v is increased within the plateau regime, localized 
states disappear one by one from inside the disordered 
domain, and those slightly above the mobility edge are 
seen to lie around its periphery. When the disordered 
domain gets well depopulated, i.e., for 40 < N v < 50, 
jx(x,y) almost vanishes inside it and there emerges a 
current circulating along its periphery, as seen clearly 
around y ~ 3 £ and y ~ 17 £ in the N v = 45 case of Fig. 
(b). This current is a diamagnetic current associated 
with the extended states surviving in the regions < y < 
5 £ and 16 £ < y < 22.5 £ of the sample, clearly seen from 
the density profile p(x, y) = ^ Q p^ a '{x, y) at x = 16 £. 

In the plateau regime 4 < N v < 48, the Hall cur- 
rent remains almost constant in net amount but changes 
in distribution with N v . When the disordered domain 
becomes well depopulated, jHaii(a;,y) virtually vanishes 
inside it and flows along its periphery in the same di- 
rection, forming two localized distributions that become 
prominent as A^ v — ► 48. See the current distribution for 
N v — 45 and its profile at x = 16 £, shown in Fig. |3|. Here 
we see explicitly that the Hall current is expelled out of 
the disordered portion of the sample. We have verified 
that even for smaller impurity strengths = 0.1 and 
Si = 0.01 the current distributions show essentially the 
same characteristics. 



V. RANDOMLY-DISTRIBUTED IMPURITIES 

In this section we examine Hall samples with ran- 
dom impurities. One of the samples we have consid- 
ered has length L x = \/2~k£ x 28 ~ 70£ and width 
L y = \/2n£ x 8 w 20£, supporting 224 electron 
states in the n = subband with < yo < 20£. 
180 short-range impurities of varying strength within the 
range |s,-| < 0.1 are randomly distributed over the full 
length and part of the width, < y < 15£. The re- 
maining region 15£ < y < 20£ is left impurity- free as 
before. We have employed a number of impurity config- 
urations for this sample and obtained qualitatively the 
same conclusion. Here we select one configuration which 
is generated by reshuffling of a basic configuration of 45 
impurities and which therefore is readily written out; see 
the appendix for details. 

Let us formally identify an electron state as localized 
if it carries a calculated amount of Hall current less than 
5 % of the typical unit (-eSE/B)(l/L x ). Then, for this 
sample 55 localized states are identified, with 10 states of 
positive energy lying in the region 9.9£ < y cm < 14.2 £ of 
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the sample width and 45 states of negative energy lying 
in the region 3.2 £ < y cm < 14.9 £. [Remember that in our 
energy scale "negative-energy" states are those that are 
captured by attractive impurities.] This contrast in both 
number and distribution of positive- and negative-energy 
localized states is direct evidence that the sample edge 
and disorder combine to efficiently delocalize electrons, 
especially positive-energy ones. These localized states lie 
within the energy range —0.128 a; < e < 0.147a;, above 
which lie only 13 edge states with y cm < 1.63 £. 
The spread in y Q of each eigenmode a 

A Q = ^|« fc | 2 |y [A : ]-(y ) a |, (5.1) 
k 

with (yo) a = X^fc l a fe| 2 ?/o[^]j is a good measure to see the 
influence of disorder on each mode. In Fig. |i| (a) we plot 
A a as a function of y cm for each state. It is clear that the 
edge states are readily identified by their small spread. 

Figure ^ (b) shows the total Hall current Juan plotted 
as a function of N v . There the formation of the upper 
plateau for N v < 90 is incomplete because of a relatively 
small number of positive-energy localized states. A steep 
fall of JHaii around N v ~ 100 implies that a large portion 
of the Hall current is carried by the 98th^ 100th states 
residing in the sample bulk 8 £ < y cm <■ 14 £, each carry- 
ing 10 ~ 40 units of (-eSE/B)(l/L x ). The subsequent 
linear decrease of Juan for 103 < iV v < 153 is due to grad- 
ual evacuation of extended states from the impurity-free 
region 15 £ ;$ y < 20 £. Finally the negative-energy local- 
ized states give rise to a clear lower plateau for iV v > 154. 

To examine the current distributions one has to recall 
that weaker impurities, as long as they trap electrons, dis- 
turb jnaii(x,y) more violently. Correspondingly, direct 
inspection of jnnii(x, y) does not clearly reveal the global 
flow of electrons. Fortunately we find the ^-averaged cur- 
rent 

JHaii(y) = j dxm^{x,y) (5.2) 

particularly suited for this purpose. 

Figure || shows the x— averaged current j™(y) — 
(1/L X ) J dxj x (x,y) and density p av {y) across the sample 
width. The paramagnetic component carried by the edge 
states is prominent in j^ v (y) only for < N v < 5. For 
N v > 10, jx V (y) is dominated by the orbital diamagnetic 
components that flow in opposite directions at opposite 
edges. The current distribution formed around y ~ 16 f is 
also a diamagnetic current associated with the extended 
states surviving in the region 15 £ < y < 20 £, seen explic- 
itly from p av (y). The density distribution p av (y) reveals 
an important role of disorder: Impurities work to dis- 
tribute electrons rather evenly throughout the sample; 
in their absence the region < y < 51, e.g., would be 
completely evacuated for N v > 70. 

The Hall current (v) scarcely deviates from a uni- 
form distribution (of the impurity- free case), drawn with 



a thin real curve in Fig. ||, while only edge states be- 
come vacant with increasing N y , i.e., for < iV v < 13. 
For N v > 14 the current distribution changes with N v , 
especially in the sample bulk. In the plateau regime 
JhbIiG/) t en ds to diminish on average in the sample bulk 
5 £ ^ y ^5 15 £ while forming prominent concentrations 
near the sample edges; see Fig. 

For comparison it is enlightening to look into the case 
where all the impurities are made repulsive, Sj — > |s,|. In 
this case no negative-energy states arise and only 3 states 
with y cm > 10 £ remain localized. No Hall plateaus are 
expected a priori. However, as shown in Fig. ^| (a), Jnaii 
varies somewhat irregularly with N v and a vague plateau 
emerges as an envelope. This implies that, while there 
are only few well-localized states, a sizable fraction of 
electron states is nearly localized, each carrying a reduced 
amount of Hall current. As in the original case, there is 
again a range of N v where juLn(y) gets small on average 
in the sample bulk; see Fig. [7] (b). 

Let us go back to the original case again and try to 
vary the Hall field. In Fig. || (a) we show the Juan — N v 
characteristics when £L is made 100 times stronger. A 
comparison with Fig. S (b) reveals an apparent decrease 
(~ 15) in the number of the electrons forming the lower 
plateau; this demonstrates that the Hall field works to 
delocalize electrons in the sample bulk. These liberated 
electrons carry a Hall current so that the steep fall of 
^Haii around N v ~ 100 becomes somewhat tempered. 
Still no drastic changes arise in the current and density 
distributions. 

However, when E y is made 1000 times stronger, i.e., 
—e£E y /uj = 1/100, delocalization prevails and only 8 
states of negative energy with y cm > 8.6 £ remain local- 
ized. An approximate linear variation of Jijaii with 7V V 
in Fig. H (b) indicates that most of the states indeed stay 
extended. Still the currents are disturbed and there is 
a general tendency that i^ji(y) concentrates near the 
two edges, though leaving a considerable amount in the 
sample bulk as well; see Fig. || (c) . 

Similarly, when E y is made 1000 times stronger, the 
vague plateau in Fig. |?] (a) disappears and a sizable 
amount of Hall current is seen to flow in the sample bulk. 
It is now clear that the decrease of the Hall current in the 
sample bulk is correlated with the dominance of localized 
states there. Figures ||, 0(b) and ||(c) thus give evidence 
that in the plateau regime the Hall current flows with a 
tendency to avoid the disordered sample bulk. 

To further elucidate the effect of the sample edges let us 
now widen the sample width to L y — ^/2tt£ x 10 w 25£, 
and embed the present configuration of 180 impurities 
within the region 5£ < y < 20 £ of the sample, leaving 
the two edge regions < y < 5£ and 20 £ < y < 25 £ 
impurity-free. 

Of the 280 electron states accommodated in the n = 
subband, we find 146 localized states lying within the 
same energy range — 0.128lj < e < 0.147a; as be- 
fore, with 73 positive-energy states lying in the region 
6.5 1 < y cm < 19.9 £ of the sample width and 73 negative- 
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energy states lying in the region 5.7 1 < y cnl < 19.8 1. 
Introduction of 56 electrons in the impurity-free buffer 
zone of width ~ 5 1 has thus brought about 91 extra lo- 
calized states. This demonstrates again how efficiently 
the sharp edge and disorder combine to delocalize elec- 
trons near the sample edges. 

Since most of the electrons are localized in the the 
sample bulk, there emerge clear Hall plateaus, as shown 
in Fig. ^ (a); there the upper plateau is not completely 
flat simply because the edge states carry a small amount 
of Hall current. Also emerges over a wide range of iV v 
the tendency that the Hall current is expelled out of the 
disordered region; see Fig. || (b); in general, as N v ap- 
proaches the mobility edge N v ~ 140, the variation of 
.?Hail(y) increases in magnitude. 

In Fig. H (c) we plot the net amount of Hall current 

per state, Jjj*J\v as a mnc ti° n °f J/ cm f° r each state. It is 
clearly seen that electron states residing on the edges of 
the disordered region ( "bulk edges" ) support a consider- 
able amount of Hall current per state and that a small 
number of states in the inner bulk carry an even larger 
amount. 

The density distribution p av (y), shown in Fig. |l^ (a), 
also reveals that in the upper-plateau regime (N v ~ 100) 
more electrons survive on the bulk edges than in the sam- 
ple interior. This feature scarcely changes even if E y is 
made 100 times stronger. When E y is made 1000 times 
stronger, impurities lose their ability to redistribute the 
electrons, as seen from Fig. [Tc] (b), and the Hall plateaus 
disappear. 

Whether localization dominates in the sample bulk or 
not depends on the competition between disorder and 
the Hall field. For the cases we have discussed so far Hall 
plateaus disappear when \si / (e£ E y / uj)\ < O(10); note 
the estimate fl3.2| ) of energy cost. For an independent 
check we have made all the impurities 1000 times weaker 
Si -> Si/1000 - O(l/10 4 ) with ~e£E v /uj = 1/10 5 kept 
fixed, and verified that the plateaus in Figs. || (b), [?] (a) 
and ^| (a) disappear with the current distributions under- 
going characteristic changes. This confirms again that it 
is the localization of electron states in the sample bulk 
that is essential for the quantum Hall effect. 

VI. SUMMARY AND DISCUSSION 

In this paper we have studied current distributions 
for finite Hall samples with disorder. The observations 
drawn from our numerical experiment are summarized as 
follows: 

(i) The current j x and the Hall current jHaii are sub- 
stantially different in distribution, (ii) The fast-traveling 
edge states carry a large amount of current per state but 
support only a negligible portion of the Hall current. In 
the Hall-plateau regime the edge states (of the uppermost 
level) are vacant and scarcely affect the distributions of 
both j x and jHaii- The edge states are in no sense the 



principal carriers of the Hall current, (iii) The sample 
edge and impurities combine to efficiently delocalize elec- 
trons near the edge, especially those to be trapped by 
repulsive impurities. Such extended "bulk-edge" states 
arise over a scale larger than the width (~ 3£) of the 
edge-state channels, (iv) In the plateau regime the Hall 
current jHaii tends to diminish on average in the sample 
bulk and concentrates near the sample edges, (v) In the 
transient regime between Hall plateaus the major por- 
tion of the Hall current is carried by a relatively small 
number of electron states in the sample bulk, (vi) The 
Hall field competes with disorder to delocalize electrons 
in the sample bulk. We have seen, in particular, that the 
breakdown of the quantum Hall effect is caused by de- 
localization of electron bulk states beyond some critical 
Hall field. 

Observations (iii) and (iv) give evidence for the expul- 
sion of the Hall current out of the_disordered sample bulk, 
which is expected theoreticallyEJ on the basis of exact 
compensation of the Hall current among extended and lo- 
calized states. This "bulk-edge" Hall current is very clear 
for the case of an array of impurities in Sec. Ill; such an 
impurity array may be realized by use of quantum 4Qts-i 
Large edge-channel widths reported experimentallytjtS 
appear to be in favor of the bulk-edge Hall current rather 
than the current carried by the edge states. 

For an improved analysis the uniform field E y we have 
employed may be replaced by a self-consistent Hall field 
generated internally by an injected current. This would 
still leave all our observations (i)^(vi) essentially intact, 
because the key factor there is competition of the Hall 
field with disorder rather than its distribution; in this 
sense the bulk-edge channels of our concern differ from 
the edge channels generated by redistributed charges £3 
Actually, the Hall field has been observed^ to be stronger 
near the sample edges than in the interior, in conformity 
with some of theoretical calculations £3 In view of obser- 
vation (vi), adopting such a realistic Hall-potential dis- 
tribution will thus work to make the bulk-edge channels 
even wider than in the uniform-field case. 

Observation (vi) suggests that field-induced derealiza- 
tion of electrons in the sample bulk is a possible mecha- 
nism for the breakdown of the QHE. With typical values 
lu ~ 10 meV and £ ~ 100 A, we find a critical field of 
magnitude \E"\ ~ 100 V/cm for the samples examined; 
in addition, keeping impurity strengths fixed yields 
the power-law dependence of Ey 1 on the magnetic field, 

E" cx £~ 3 cx £? 3 / 2 . A critical field of this order of magni- 
tude and the B dependence of this form appear consistent 
with experiments. E3 It would be worthwhile to further 
explore this mechanism. 
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APPENDIX A: IMPURITY DISTRIBUTION 



In this appendix we display the impurity configuration 
used in Sec. V. Let us first select 45 points {(x~i,y~i);i — 
1, • • • , 45} randomly distributed over the range < x < 6 
and < y < 6. One such choice we adopt is 
P 45 ={ (0.28, 3.85), (0.34, 0.33), (0 
(0.57, 5.52), (0.64, 1.23), (0.79 



(0.95, 
(1.68, 
(2.20, 
(2.77, 
(3.34, 
(3.72, 
(4.55, 
(4.85, 
(5.35, 



0.61), 
3.66), 
1.62), 
1.14), 
4.01), 
0.26), 
1.58), 
2.88), 
3.80), 



(1.05, 
(1.77, 
(2.47, 
(2.81, 
(3.41, 
(3.91, 
(4.61, 
(4.90, 
(5.46, 



56, 2.32), (0.40, 5.04) 



4.78), (0.91, 5.86) 
3.70), (1.22, 2.75), (1.41, 4.43) 
0.73), (1.86, 0.24), (2.04, 5.32) 
2.29), (2.55, 5.46), (2.72, 2.86) 
4.46), (2.93, 3.70), (3.33, 1.83) 
2.83), (3.64, 5.15), (3.68, 0.99) 
3.40), (4.33, 4.75), (4.37, 4.19) 
3.32), (4.65, 0.64), (4.76, 4.12) 
5.08), (5.03, 5.48), (5.21, 1.57) 
1.11), (5.57, 2.42), (5.60, 4.66) 
(5.96, 0.16) }. 

Construct out of this P 45 another table {(y~i + 6, xi); i = 
1, • • • , 45} (via interchange Xi <-> y% and a shift), and sort 
it into ascending order with respect to the + 6 compo- 
nents to form P\ 5 = {(6.16, 5.96), (11.86, 0.91)}. 
These P45 and P45 are combined to yield a table of 
90 points, P 90 = {(0.28, 3.85), (11.86, 0.91)}, dis- 
tributed over the range < x < 12 and < y < 6. 

On the other hand, we shift P45 by (3,3) (which 
is our arbitrary choice) modulo 6 and reorder the re- 
sulting table {(xi + 3,y~i + 3) mod 6; i = l,---,45} 
into ascending order with respect to the first compo- 
nents to form P'45 = {(0.33, 4.83), (5.93, 0.7)}. 
We then construct out of this P'45 a new table P'90 = 
{(0.33, 4.83), • • • , (11.88, 1.85)} of 90 points by the same 
procedure as done for Pg out of P45 . 

Finally wc shift P' 90 by (12,0), (£»,&) -> {xi + 12, yi), 
and add it to P90, obtaining the table of 180 points 
Piso = {(0.28, 3.85), •••,(23.86, 5.72), (23.88, 1.85)} 
distributed over the range < x < 24 and < y < 6. Via 
the rescaling Xi = V2n£ x (28/24) x Xi and = ^/2tt£ x y~i 
we distribute these 180 impurities at positions (xi,yi) 
within the region < x < 70£ and < y < 15 £ of the 
sample. 

The strengths s t we adopt for the first 45 impurities 
are 

S 45 = O.lx {-0.973, -0.826, 0.855, -0.368, 0.793, 0.316, 
-0.233, 0.180, -0.692, 0.916, -0.463, 0.252, -0.452, -0.745, 
-0.485, 0.446, -0.622, -0.076, 0.503, 0.205, -0.973, 0.758, 
0.544, 0.720, 0.997, 0.684, 0.689, 0.088, -0.874, -0.631, - 
0.078, 0.908, 0.818, -0.548, -0.6141, -0.344, 0.271, -0.802, 
0.871, 0.218, -0.106, 0.283, -0.632, -0.986, -0.132}, 
which are generated randomly within the range \ si\ < 0.1. 
To generate 180 strengths we arbitrarily choose to set 
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FIG. 1. Electron states localized about isolated impu- 
rities and an extended state. (a) Density distributions 
p^ a \x,y) plotted in units of 1/(£L X ), (b) Current distribu- 
tions jx {x,y) in units of —eu/L x . (c) Hall-current distribu- 
tions jy"Ju(x,y) in units of — (e 2 /2irTi)5E y . The coordinates 
x and y are plotted in units of the magnetic length I. 




FIG. 2. (a) Currents vs N v . The total Hall current Jhm 
is plotted in units of —(eSE y /B)(l/L x ), and J x in units of 
0.05 x (— euil/ L x ). (b) Current and density profiles at x — 161 
for iVv = 0, 10, 45, with j x (x, y) plotted in units of —eLu/(2nl) 
and p{x,y) in units of l/(27r^ 2 ). 
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FIG. 3. (a) The Hall current density jna,n(x,y) for the 
n = level with N v = 45, plotted in units of —(e 2 /2ivh) 5E y . 
(b) Current profile at x = 
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FIG. 4. A sample with random impurities, (a) The spread 
A a of the wave function in yo space (in units of £) plotted 
as a function of y cm for each state, (b) Currents vs iV v . 
The JHaii and J x are plotted in units of —(eSE y /B)(l/L x ) 
and 0.05 x (~ecu£/L x ), respectively, with 8E y /E y = 0.01 and 
e£Ey/u = -1/10 5 . 



0.2; f p~N v =W 

o.iA'. ; ' v -r^w v =7o 



.av 

J x (y) 




-N y =l00 



-1-0" 



15 




P (y) 



-jv v =ioo 



5 10 15 20 y 

FIG. 5. Current and density distributions for the n — 
subband, with j x v (y) plotted in units of —eu/(2nl). 




FIG. 6. Hall-current distributions Jh^u(j/), in units of 
-(e 2 /2-kTi) SE y , for the n — subband in the plateau regime. 
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FIG. 7. (a) Currents vs iV v and (b) Hall-current distribu- 
tions in the case of random repulsive impurities, Si — ► \si\. 
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FIG. 8. (a) Hall plateaus persist when E y is made 100 times 
stronger, (b) Plateaus disappear when E y is made 1000 times 
stronger, (c) A considerable amount of Hall current thereby 
flows in the sample bulk. 
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FIG. 9. A Hall sample with an impurity-free buffer zone of 
width ~ 5£ along the y — edge, (a) Currents vs N v . (b) 
Hall-current distributions, (c) Net amount of Hall current per 
state [on a square-root scale in units of — {e5E y / B)(\ / L x )\ as 
a function of y cm ; here omitted for ease of exposition are 3 
pairs of points (mostly out of the range shown) that combine 
to carry only -7.8 units of current. 
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FIG. 10. Density profiles for the n — subband with 
iVv=0,50,100,140 and 160. (a) e£E y /io = -1/10 5 and (b) 
elEyjuj = -1/100. 
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